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We calculate the intersubband absorption linewidth 2rop in quantum wells (QWs) due to scat- 
tering by interface roughness, LO phonons, LA phonons, alloy disorder, and ionized impurities, and 
compare it with the transport energy broadening 2rtr = 2ft/rtr, which corresponds to the transport 
relaxation time rtr related to the electron mobility ^. Numerical calculations for GaAs QWs clarify 
the different contributions of each individual scattering mechanism to the absorption linewidth 2rop 
and transport broadening 2rtr. 

Interface roughness scattering contributes about an order of magnitude more to the linewidth 2rop 
than to the transport broadening 2rtr, because the contribution from the intrasubband scattering 
in the first excited subband is much larger than that in the ground subband. On the other hand, LO 
phonon scattering (at room temperature) and ionized impurity scattering contribute much less to 
the linewidth 2rop than to the transport broadening 2rtr- LA phonon scattering makes comparable 
contributions to the linewidth 2rop and transport broadening 2rtr, and so does alloy disorder 
scattering. 

The combination of these contributions with significantly different characteristics makes the ab- 
solute values of the linewidth 2rop and transport broadening 2rtr very different, and leads to the 
apparent lack of correlation between them when a parameter, such as temperature or alloy compo- 
sition, is changed. Our numerical calculations can quantitatively explain the previously reported 
experimental results. 

PACS numbers: 78.67.De, 78.30.Fs, 73.21.Fg, 73.63.Hs. 



I. INTRODUCTION 



The intersubband absorption linewidth in semiconductor quantum wells (QWs) closely relates to fundamental 
problems in the physics of optical transition, such as relaxation many-body effects |j, and disorder |j. 
Furthermore, it is a key factor in improving the performance of quantum cascade lasers |^ and QW infrared photode- 
tectors 

To investigate the effects of various scattering processes, intersubband absorption linewidths have been measured 
for various temperatures well widths alloy compositions and doping positions in GaAs and other QWs. 
These results show that absorption linewidth has a weak dependence on temperature and alloy composition and 
apparently has little correlation with mobility. Its strong well-width dependence suggests that the main contribution 
is from interface roughness scattering. 

In a previous paper [ pd] , we discussed the effect of interface roughness scattering on linewidth by comparing 
calculations based on a microscopic theory by Ando and experimental data for modulation-doped GaAs/AlAs 
QWs with a well width of 80 A. The results made it clear that linewidth is much more sensitive to interface roughness 
scattering than transport mobility is, because the contribution from the intrasubband scattering in the first excited 
subband is much larger than that in the ground subband Even in wide GaAs QWs, where interface roughness 
scattering should be less effective, recent reports [12, 1^ showed that interface roughness scattering has a larger effect 
on linewidth than either electron-electron scattering or bulk impurity scattering. 
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In the present paper, we apply our theoretical method 1 11 to scattering by LO phonons, LA phonons, alloy disorder, 



and ionized impurities as well as interface roughness scattering, in order to compare their respective contributions 
to intersubband absorption linewidth and transport mobility. Numerical calculations for GaAs QWs confirm that 
the very high sensitivity of linewidth to interface roughness scattering is the key to quantitatively explaining the 
previously reported experimental results for linewidth in comparison with mobility. 

The method presented here follows Ando's theory [||, in which the intrasubband and intersubband energy- dependent 
single-particle [ [4| relaxation rates for various scattering mechanisms are first calculated and then included in a formula 
for the two-dimensional (2D) dynamical conductivity Re fizz (w) to give the absorption lineshape for photon frequency 
LO. This method is similar to a familiar method of calculating transport mobility |15, |l^. It is important to note 
that intersubband optical absorption is the collective excitation among a confined electron gas. However, our present 
calculation of single-particle relaxation rates and lineshape is very important and useful, because the absorption 
lineshape Re (122(0;) of collective excitation is given by the single-particle dynamical conductivity azz{^) via p7 

) 

-C^zzi^) 



azz{uj) = (1) 



in the crudest approximation. Here, is the vacuum permittivity, kq is the static dielectric constant of the 2D 
material, and dcff is the effective thickness of the 2D electron gas [|l5[ . 

The collective excitation effects, or many-body interaction effects, on intersubband absorption linewidth have been 
issues of recent interest in both theoretical and experimental studies. In the limit of small band-nonparabolicity and 
constant single-particle relaxation rates, Nikonov et al. theoretically showed that many-body effects only cause blue- 
shifts in absorption spectra (the depolarization shift) and that the linewidth is solely determined by the single-particle 
relaxation rate |Q . 

In largely nonparabolic systems, the variation in energy separation between the ground and first excited subbands 
produces additional width in the single-particle excitation lineshape RetTzz(a;). However, many-body effects lead to 
redistribution of oscillator strength and collective excitation that has a sharp resonance. As a result, the linewidth 
of the collective excitation spectrum Rea'zz(aj) is significantly different from that of the single-particle excitation 
spectrum Re i7zz(ijj) |l^, Furthermore, nonparabolicity causes difficulties in calculating single-particle relaxation 
rates and Re 1722(0;) by Ando's formalism |20|. Experiments to elucidate these effects were performed by Warburton 
et al. on InAs/AlSb QWs [|. 

In the more popular systems of GaAs/AlGaAs and InGaAs/InAlAs QWs, however, our present calculation, which 
assumes small nonparabolicity, is applicable. The purposes of this paper are (1) to calculate intrasubband and 
intersubband single-particle relaxation rates for relevant scattering mechanisms as functions of in-plane kinetic energy 
assuming small nonparabolicity and (2) to quantitatively explain previously reported experimental data on linewidth 
and mobility in GaAs-based QWs, which appeared to have little correlation. 

In the next section, we summarize how linewidth and mobility are related to single-particle relaxation rates, and 
calculate the single-particle relaxation rates for various scattering mechanisms as functions of the kinetic energy E. It 
is shown that linewidth and mobility have very different sensitivities to the same scattering mechanism. In Section 3, 
previously reported experimental data for various temperatures, well widths, alloy compositions, and doping positions 
are quantitatively explained by numerical calculations. 

II. FORMULATION OF THE PROBLEM 
A. General theory of intersubband absorption lineshape and transport mobility in quantum wells 

A general theory of intersubband absorption linewidth due to elastic scatterers in 2D systems was formulated by 
Ando ^ |2^. According to Ando's theory, the absorption lineshape of single-particle excitation between the two 
lowest subbands can be expressed as 

Reazziuj) = / —^dEf{E)— ^2 .p^2 ' 2 

when all electrons are initially in the ground subband. Here, 

^op{E) = i[ri„tra(£^)+ri„tcr(£^)], (3) 

T,,,UE) - 2^^(|(0k'|ifi|0k)-(lk'|i/i|lk)nj(e(k)-e(k'))|^^^(^), (4) 



3 



r^nUE) = 2^^(|(Ok'|ffi|lk)n5(e(k)-£(k')+£;io)U^,(k)' (5) 

k' 

e is the elementary charge, h is the reduced Planck constant, m* is the electron effective mass, /lo is the oscillator 
strength, Eio (= Ei—Eq) is the intersubband energy separation, f{E) is the Fermi distribution function at temperature 
T, |nk) is the state vector of the electron with subband index n and wave vector k. En is the quantization energy, 
e(k) = fi^k'^ /2m* , Hi is the scattering potential, and (• • •) denotes the average over distribution of scatterers. This 
theory assumes a parabolic conduction band, or a constant effective mass for different subbands; a modification for 
slightly nonparabolic systems like GaAs QWs will be described in a later paragraph. In this paper we denote the full 
width at half maximum of the spectrum given by Eq. as 2rop. 

Note, on the other hand, that the transport relaxation time Ttx{E), or the transport relaxation rate 2rtr(-E) — 
'2-Ti/Ttt{E) can be expressed as 

= 47r5] ( |(0k'|i7i|0k)p ) 5 (e(k) - £(k')) (1 ~ cosB) |^^^(^^ , (6) 

k' 

where is the angle between k and k'. The mobility is given by // — er^-i-lm* with an average relaxation time of 

--/--(-)-^//--^- 

To enable quantitative comparison between the linewidth 2rop and mobility /i, we define the transport energy 
broadening as 2rtr — 2h/Ttr — 2fie/m*^ In particular, low-temperature transport broadening is given by 

2rtr = 2Ttr{EF) = 2h/Tty:{EF), whcrc Ep is the 2D Fermi energy. 

There are two relevant many-body effects: static and dynamic screening. The former screens the potentials of 
elastic scatterers while the latter induces collective charge-density excitation because of the incident optoelectric field. 

The static screening effect can be included by replacing the scattering matrix element (mk'|iJi|rik) with W 



(mk'|i?i|nk) + (0k'|i7i|0k) 



1 

- 1 



F(00)(mn){q) , . 

-f'(oo)(oo)(,9j 



Here, q = k — k', e{q,T) is the static dielectric function ||T^, |T^, and i^(fcz)(mn) is a form factor defined by ||l| 

Fiki)imn){q) ^ I dz ( dz'a(^)0(^)C™(^')Cn(^')e~«l^-"'l. (9) 



The z axis is set along the growth direction of samples, and Cn(^) is the wave function for the n-th subband electron 
motion in the z direction, which is chosen to be real. The screening correction only results in dividing (Ok'|iJi|Ok) in 
Eq. (I) by the factor 



1 ( ^ _] \ Fm(ii){q) 



<^{<l,T) \e{q,T) J i^(oo)(oo) (?) 



(10) 



and (Ok'|_ffi|Ok) in Eq. (^ by e{q,T). In this paper we only treat symmetrical QWs, so there is no screening factor in 
Eq. (H). e{q, T) significantly increases mobility, while ^(g, T) hardly affects absorption linewidth because S{q^ T) ^ 1. 

The dynamic screening effect is counted as a depolarization field, and the absorption lincshape Re(Tzz((jj) of the 
induced collective charge-density excitation is given by 



crzz{(^) = ^ (11) 



with the dynamical dielectric function of 



e2,(cj) = l + —azzicj). (12) 

The resonance energy Eiq of Kcazzii-^) is blue-shifted from the original resonance energy Eiq of KecTzzi'-^), and 



Eio = V^io' + i^^p)^ (13) 
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with the plasma frequency of 



V eo^o™ "cff 

The blue-shift Eiq — Eiq « (Tiujp)'^ / (2£'io) is called the depolarization shift. The linewidth of Re Hzz {(^) is the same as 
that of Rc azz (w) if 2rop(i?) is independent of energy ^ , though they are different in general. When the depolarization 
shift is sufficiently small, or 

^10 - £^10 < 2rop(0), (15) 

(Jzzi'-^) is approximately equal to ^^^(w). 

Although Eqs. (|)-(|) were derived assuming parabolic bands, we may apply them to slightly nonparabolic systems 
in which the additional width due to nonparabolicity is small compared with the width due to scattering mechanisms. 
The condition is expressed as {l—mQ/mDEp < 2rop(0) at low temperatures, where m* is the electron effective mass in 
the n-th subband. In this case, we can use the present theory by replacing Eiq in Eq. (^ with £^io(0) — (1 — mQ/ml)E 
p9| , [2l} |, which has a much larger influence on absorption linewidth than other corrections. Here, i?io(0) represents 
the intersubband energy separation at k = 0. For consistency, respective i5-functions appearing with the squares 
of scattering matrix elements |(0k'|ffi|0k)|2, |(lk'|ifi |lk)|2, |(Ok'|Hi|Ok)(lk'|ffi|lk)|, and |(0k'|iJi|lk)|2 in Eqs. 
(|) and (|) should be replaced by (5(eo(k) - eo(k')), S{ei(k) - ei(k')), ^ [(5(eo(k) - £o(k')) + S{ei{k) - £i(k'))], and 
(5(£i(k) — £o(k')), where e„(k) = En + ?i^fc^/2m*. Values of m* obtained from the Kane model are used in our 
numerical calculations. 

In most cases of GaAs QWs examined later in this paper, the depolarization shift Eiq — Eiq and the nonparabolicity 
effect (1 — mQ/mVjEp are small compared with 2rop(0), so the absorption linewidth is estimated directly from 
Rea^,(c^) inEq. @. 



B. Scattering mechanisms 

In this section, we calculate and compare 2rop(£') and 2rti.(i?) due to scattering by interface roughness (IFR), 
LO phonons, LA phonons, alloy disorder (AD), and ionized impurities (ION). Furthermore, numerical calculations of 
each individual scattering mechanism are performed for modulation-doped GaAs (or InGaAs)/AlAs QWs. In actual 
samples, several scattering mechanisms coexist; the total scattering rate can be obtained as the sum of their rates. 
Namely, 

T,^{E) = r(ip™) [E) + r(LO) (^E) + r(LA) {^e) + F'^d) [E) + V^^''^ {E) + ---, (16) 

r,AE) = F|r^(i?) + F|^^°)(i?) + r^^^\E) + Fr^i?) + rii°^^(i?) + • • • . (i7) 

For simplicity, we perform numerical calculations in single QWs with a finite barrier height of Vq where band 
bending due to doping is neglected. The origin of the z axis is set at the center of the QWs. Material constants of 
GaAs used in calculations are shown in Table 0. 



1. Interface roughness scattering 

In GaAs/AlGaAs QWs, dominant monolayer (ML) fluctuations are formed at the GaAs-on-AlGaAs interface (Al- 
GaAs surface covered by GaAs). We assume that the roughness height A(r) at the in-plane position r = {x,y) has a 
correlation function [ITsl : 



(A(r)A(r')) = cxp [-^^^) , (18) 
where A is the mean height of roughness and A is the correlation length. The scattering matrix element is given by 

(mk'|Hi|nk)= J d\F^n ^{r) e'^^ " (19) 

with 

F„^n = Vo Cmi-L/2) Cn(-i/2), (20) 
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where L is the well width and C„(— L/2) is the wave function at the GaAs-on-AlGaAs interface. Because interface 
roughness is equivalent to local fluctuations in well width, i^mn in Eq. ( |20| ) can also be expressed as 



Fr,^n = VidE„JdL){dE„/dL). 
In the case of the infinite-barrier approximation, Eq. ( pO| ) can be expressed in an alternative form as [ p5| p3| 

h'^ dCrn{z) dC„{z) 



p — 



2m* dz dz 



(21) 
(22) 



which is found to be proportional to L ^ . 

Substituting Eq. (|l|) into Eqs. (|) and (|), we get 



^ intra \^ I 



2 a2 /-Tr 



.(IFR) 



h JO 
m*A2A2 



de 



z=-L/2 



Fu 



S{q,T) 

dee-^'^'i\ 



where the absolute values of the 2D scattering vectors q and q are given by |l] 



2fc2(l 



2k^ 



2m* El 



2k\ k 



2m* El 



On the other hand, we can express the transport relaxation time t^I^^'\e) |2^, or the transport relaxation rate 



■ cos 9. 



(23) 
(24) 

(25) 
(26) 



2rf«)(i?) = 2n/rr"^(£;) as 



(IFR) , 



2r. 



(IFR) 



{E) = 



2m*A2A2 



de- 



1 — cos ( 



(27) 



which is similar to Eqs. (|2^) and (^J). 

It is useful here to comment on the similarities and differences in the equations for rintra(£^), rinter(£^), and 2Ttr{E). 
First, all three are proportional to A^, and also to A^ for small A. Second, Tint^y{E) is much smaller than rintra(-£'), 
because q is smaller than q. Third, {1 — cos 9) / e{q,T)'^ appearing in 2rtr(£') shows that the forward scattering {9 ~ 0) 
does not contribute to transport broadening, and that the screening effect reduces the scattering rates. Finally, and 
most importantly, they include different factors [Foo/S{q,T) — i^n]^, Fqi"^, and 2Foo^. S{q,T) can be neglected 



because Ei is more sensitive to L than Eq- 
) As a result, rintra(-E') is much larger than 



because S{q,T) ^ 1. As is shown below, Fu is much larger than Fqq, 
(In the infinite-barrier approximation, Fu is four times larger than Fqq 

2ruiE). 

Figure shows 2rop(£'), rintra(-B), rintcr(£'), and 2rtr(-B) in a modulation-doped GaAs/AlAs QW with L = 80 A, 
A = 3 A, and A = 50 A. These values of A and A are typical for the GaAs-on-AlAs interface ||l^, |2^ . Temperature 
was set at T = OK, and sheet electron concentration was chosen to be Ns = 5 x 10^^ cm~^, which gives Fermi energy 
of Ep — 17.8 meV. The same values of L, Ns, A, and A are also used for calculations of other scattering mechanisms 
in this section. 

In Fig. 1^, rintra(-£') dccrcascs as E increases, and it has a maximum value of 8.3 meV at E — OmeV. rintcr(-E') 
is almost constant with respect to E, and its value of 0.6 meV is much smaller than that of rintra(^')- The values 
of 2rtr(F) are OmeV a,t E = OmeV owing to the screening effect, and 0.6 meV at E = Ep, which determines the 
low-temperature transport broadening. As a result, 2rop(i?), the sum of rintra(-E) and rintcr(-£'), is found to be much 
larger than 2rtr(£'). 

Figure || shows 2rop and 2rtr as functions of correlation length A; they are calculated respectively from Eqs. (|^) 
and (1) in a modulation-doped GaAs/AlAs QW with L = 80 A, iV^ = 5 x 10" cm-2, T = K, and A = 3 A. 2rop,para 
represents the linewidth calculated without changing Eio into ^10(0) — (1 — m^/ml)E in Eq. (||). 



In Fig. H, 2ro 



and 2rtr are both proportional to A for small A with the difference in absolute values being 



about one order of magnitude. With nonparabolicity, 2rop has a lower limit of {l — mQ/mDEp = 1.35 meV in addition 
to 2rop,para. For large A, the insensitivity of 2rtr to forward scattering causes its value to be smaller. This shows 
that the correlation length of A ^ l/kp contributes most to 2rtr, where kp is the Fermi wavenumber. In principle, 
values of the roughness parameters A and A can be uniquely determined if linewidth and mobility arc both measured 
at low temperatures. 

As shown above, the main characteristic of interface roughness scattering is its order-of-magnitude different contri- 
butions to linewidth and transport broadening (and hence to mobility). This is the key point for understanding the 
apparent lack of correlation between them. 
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FIG. 1: 2rop(-E), TintTa{E), rintcr(£^), and 2TtT{E) due to interface roughness (IFR) scattering, plotted as functions of the 
in-plane kinetic energy E. 
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FIG. 2: Intersubband absorption linewidth 2rop and the transport energy broadening 2rtr due to interface roughness (IFR) 
scattering, plotted as functions of the correlation length A. The 2rop,para represents the absorption linewidth without the 
additional width due to band-nonparabolicity. 



2. LO phonon scattering 



In considering phonon scattering processes, it should be noted that phonons have approximately three-dimensional 
(3D) properties, since they are hardly confined to QWs. The z-component momentum conservation in 3D systems 
requires the scattering matrix element of 2D electrons to be calculated from |^ 

|M2Dp = 5]|A^3D/(g.)P. (28) 

Here, M2D and M3D are the 2D and 3D scattering matrix elements, respectively, and I{qz) is given by 

I{q.) - /™„(g.) = [ dzU{z)Cn{z)e"''\ (29) 



Since the method of calculating the 3D scattering matrix element is well established, the 2D scattering matrix element 
can be easily obtained. 

In polar optical phonon scattering, or simply LO phonon scattering, the 3D scattering matrix element is given by 



(30) 
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for phonon emission processes, and by 



2\ e'^huji^oNw f 1 



2eoQ^ 



Kq 



(31) 



for phonon absorption processes. Here, Q is the absolute value of the 3D scattering vector, Koo is the optical dielectric 
constant, wlo is the LO phonon frequency, and A'lo is the LO phonon occupation given by 



1 



(32) 



Since LO phonon scattering is an inelastic process, Eqs. (^ and (^ are not applicable in their original forms. However, 
by modifying the (5-functions in Eqs. (^ and (||), 



6{e{k)-e{k')) 
(5(£(k)-e(k') + ^io) 



6 (e(k) - e(k') ± nu;w) , 
5(e(k)-£(k') + Sio±?ic^Lo) 



(33) 
(34) 



such that total energy is conserved, we can estimate width of the zero-phonon band in an approximation that neglects 
phonon sidebands. Here, ± indicates phonon absorption (+) and emission (— ). Thus, we have 



r 



(LO) 



m*e LUi^O / I 1 
Aireoh \koo kq J Jo 



de 



X 



id{E- huji^o) {F(oo)(oo)(gc) - 2f(oo)(ii)(qc) + -F(ii)(ii) (^c)} 



(iVLo) 



.(LO) 



{F{QO){00){qa) - 2i^(00)(ll)(9a) + -F'(ll)(ll) (^a) } 
d0 



m*e cjlo / I 1 
Anenfi \Koo / Jo 



X e{E + Eio - huj-Lo) ^^^"^^'' f(01)(10)(gc) + ^^^°^ f(01)(10)((?a) 

where &{E) is the Heaviside step function. Absolute values of scattering vectors are given by 

2m*tjLO 





= 2P- 




= 2fc2 + 


qc' 


= 2fc2 + 


qa' 


= 2fc2 + 



h 

2m*tJLO 

n 

2m*Eio 



2fc,/p + ^!I^cos^, 



2m*WLO 



2 , 2m *Eio ^ 2m*uJLo _ /fc2 + '^m*Eio 



2m*WLO 



COS0, 



cos 6*, 



(35) 

(36) 

(37) 
(38) 
(39) 
(40) 



and the subscripts "e" and "a" represent emission and absorption of LO phonons, respectively. 
On the other hand, the transport relaxation rate can be expressed as 



2T^^''\e) 



m*e LOi^o ( 1 
27reo?i 



1 



1 



KoJ l-f{E) Jo 



e{E- hui^o) {1 - f{E " ?iwLo)} ^^^°"'"^^ f(oo)(oo) (^c) 



+ {1- fiE + hiui^o)} 



(TV, 



LO/ 



9a 



(00) (00) 



('Za) 



(41) 



in the approximation that neglects the "in-scattering term" [|l6[ . Since LO phonon frequency is high, the screening 
effect can be neglected. 

The four form factors F(oo)(oo) (?), F{QO){ii){q), ^(ii)(ii) and i^(oi)(io) (9) appearing in Eqs. (H), and ^) 
are plotted as functions of q in Fig. ||. First, i^(oo)(oo) (?)) F{ao){ii)iq)7 and -F(ii)(ii)('i') are very close, which makes 
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FIG. 3: Dependence of the form factors -F{oo){oo) (9)1 ^'(oo)(ii)('?)i (9)) and -F{oi){io) (9) on the absolute value of the 

two-dimensional scattering vector q. 
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FIG. 4: 2rop(-B), rintra(i?), rintcr(i?), and 2Vty:{E) due to LO phonon scattering, plotted as functions of the in-plane kinetic 
energy E. LO phonon energy is fiu)-Lo — 36.5 meV. 




F{oo){oo){q) - 2-F(oo)(ii)(9) + in ri„tra(-B) much smaller than i^(oo)(oo) (?) in 2TtriE). In other words, the 

difference in intrasubband scattering matrix elements for the two subbands is small in LO phonon scattering. Second, 
■F'(oi)(io)(9) in ri„tor(^^) is much smaller than F(oo)(oo) (9)- 

Figure | shows 2rop(£'), Tintia{E), TinteriE), and 2rti.(£') due to LO phonon scattering at T = 300 K. First, 
FintraC-E') is much smaller than 2Fti(-E), because the difference in intrasubband scattering matrix elements for the 
two subband is small in LO phonon scattering as already shown in Fig. ^. Second, FinterC^') is much smaller than 
2Ttr{E) owing to the small form factor of f(oi)(io) (?) and the large absolute value of scattering vector q. Third, when 
the kinetic energy E is larger than the LO phonon energy of Ej^o — 36.5 meV, intrasubband LO phonon emission 
is allowed, which makes 2Ftr(£') and Fintra(-£') larger. As a result, 2Fop(-E) is much smaller than 2Ftr(-E) at room 
temperature. 

When systems are cooled down to OK, only intersubband LO phonon spontaneous emission is allowed in the case 
of E'lo > i?LO- Therefore, Fintra(E^) and 2Tti{E) vanish, and only TinteriE) has a finite value of about ImeV. 



3. LA phonon scattering 

Acoustic phonon scattering via deformation potential coupling, or simply LA phonon scattering, is virtually elastic. 
The 3D scattering matrix element is given by ||2^ 

(|M3dP) = ^^ (42) 
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FIG. 5: 2rop(i5), rintra(i5), rintcrl-E), and 2rtr(-B) due to LA phonon scattering, plotted as functions of the in-plane kinetic 
energy E. 



for both LA phonon emission and absorption processes, where D is the deformation potential constant and q is the 
longitudinal clastic constant. Note that Eq. (j^j is independent of the scattering vector as a result of the linear 
dispersion relation of LA phonons. Therefore, we have 



.(LA) 



.(LA) 



(E) 



m*kBTD'^ 
m*kBTD^ 



de / dz 



S{q,T) 



d0 / dz[Co{z)Ci{z)Y 



(43) 
(44) 



^intcrH-^) independent of the kinetic energy E, and ^InCii-E) is also almost independent of it because S{q,T) ~ 1. 
On the other hand, the transport relaxation rate is given by [ p5| , p5[ 

2m*kBD^T 



.(LA) 



2r 



(LA) 



(E) = 



dO 



1 — cos t 



dzC,Q{zY 



(45) 



which has an energy dependence due to the screening effect. 

Note here that the z-integrals of C,q{zY, (li{zY^ and [C,i){z)(^i{z)]'^ have comparable values (3/2L, 3/2i, and 
respectively in the infinite-barrier approximation); thus rintra(£') and ri,itcr(£') are nearly equal, and 2rop(£'), the 
sum of them, is comparable with 2rtr(£'). 

Figure I shows 2T^p{E), T-.^tUE), rintcr(£^), and 2TtAE) at T = 300 K. rintra(^^) and rintcr(£^) have almost the 
same constant values of about 0.5 meV. 2rtr(-E) vanishes a,t E = OmeV owing to the screening effect, and approaches 
a constant value of about 1.5 meV as E increases. 



4- Alloy disorder scattering 



When there are alloy layers composed of A^Bi^xC, such as Ala;Gai_a;As and In^jGai-xAs, electrons are scattered 
by conduction band disorder. The scattering matrix element due to alloy disorder is given by 123, |2^ 



(|(™k'|i7i|nk)n = ^!i^^£lMi^ 



(46) 



alloy 



where a is the lattice constant and SEc is the difference in conduction band minima of crystals AC and BC (AlAs and 
GaAs in the case of A^Gai-^As). Note that Eq. ( p6| ) is independent of the scattering vector owing to the short-range 
nature of the scatterers. Therefore, we have 



^ intra l-C/j - 



.(AD) 



(E) = 



m*a'^{6Ecyx{l- x) 
m*a^{6Ec)^x{l~ x) 



de 
de 



dz 



alloy 



Co(^)^ 



CM' 



dz [Co(^)Ci(^)]'. 



(47) 
(48) 
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FIG. 6: 2rop(-E), rintra(i?), ^iTAaiiE), and 2rtr(-E) due to alloy disorder (AD) scattering, plotted as functions of the in-plane 
kinetic energy E. 



rintor(£') is independent of i?, and rintra(£') is also almost independent of it because S{q,T) 
On the other hand, the transport relaxation rate is given by [^6|, ^ 



nh' Jo e{q,T) 



alloy 



which has an energy dependence due to the screening effect. 

Since alloy disorder scattering is due to 5Ec, it can be regarded as a kind of roughness scattering. If one substitutes 
Vo = 5Ec, = a?x{l — x^-i, and = a^/27r into Eqs. (p3|), (p^, and (p7|), one can recognize that the alloy 
disorder scattering of Eqs. (^7|)-(^9|) is expressed as the sum of the "roughness scattering" rates due to the alloy layer 
at position z. 

Note here that rintia(-E'), FintorC-E'), and 2rtr(£') f^or alloy disorder scattering are similar in form to those for LA 
phonon scattering; thus 2rop(-E) is comparable with 2rtr(£'), as in LA phonon scattering. 

Figure I shows 2Fop(-B), Fi„tra(-B), ^int^E), and 2Fti.(-B) in luo.iGao.gAs/AlAs QWs (x = 0.1) at T = OK. We 
set the lattice constant and conduction band offset to approximately a — 5.66 A and 5Ec — 700 meV, respectively. 
rintia(-E') and Fintcr(-E') have almost the same constant values of about 0.2 meV. 2Fti (-E) vanishes a.t E = OmeV owing 
to the screening effect, and approaches a constant value of about 0.57 meV as E increases. 



5. Ionized impurity scattering 

When dopant donors of Si are ionized, electrons supplied to QWs suffer from scattering by the Coulomb potential 
of the donors. The scattering matrix element due to an ionized impurity at position Z is given by j^] 



(mk'|iJi|nk) = 



2eoKo9 



(50) 



Therefore, we have 



^(lON) 



(E) = 



m*e' 



o4 



dZN{Z) / de 



dZN{Z) / d0 



qj ""15(9, T) 



dz 



-q\z~Z\ 



(51) 
(52) 



where N{Z) is the 3D impurity concentration at position Z . The transport relaxation rate, on the other hand, is 
given by pi ra 



2r 



(ION) 



(E) 



* 4 

me 



27reo^Ko^^^ 



dZNiZ) / d0 



1 — cos 9 
q'e{q,Ty 



dz(^o{z) e 



2 ~q\z-Z\ 



(53) 
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Figure | shows 2rop(i^), Tintrs.{E), TinteriE), and 2Ttr{E) in a J-doped GaAs/AIAs QW with 60 A spacer layers 
{Z = 100 A) at T = OK. First, rintra(£') is much smaller than 2rtr(£'), because the difference in intrasubband 
scattering matrix elements for the two subbands is small in ionized impurity scattering. Second, Pinter (^') is much 
smaller than 2rtr(£') owing to the large absolute value of scattering vector q. As a result, 2rop(-E) is much smaller 
than 2rtr(£'). 



III. COMPARISON WITH EXPERIMENTAL RESULTS 



In this section, the absorption lincwidth 2rop and transport energy broadening 2rtr are calculated for some GaAs 
(or InGaAs)/AlAs (or AlGaAs) QWs as functions of temperature, well width, alloy composition, and donor doping 
position. These results are compared with previously reported experimental data. 



A. Temperature dependence 

Experimental measurements of the temperature dependence of absorption lincwidth are expected to clarify the 
effects of phonon scattering. We previously reported absorption linewidths in comparison with transport mobilities 
in a modulation-doped GaAs/AIAs single QW with a well width of i = 80A and a sheet electron concentration of 
Ns = 9.8 X 10^^ cm^^, at temperatures ranging from 4.5 to 300 K [0. The absorption spectrum observed at 4.5 K is 
shown in Fig. ^ The low-temperature lincwidth 2rop was 11.1 meV and the low-temperature transport broadening 
2rtr = 2he/mQ^ was 1.2 meV, which was calculated from the mobility n of 2.9 x lO'^cm^/Vs. Note that lincwidth 
was one order of magnitude larger than transport broadening at low temperatures. The temperature dependences of 
lincwidth and transport broadening are plotted in Fig. ^ by solid and open circles, respectively. 

We performed the corresponding calculations of lincwidth and transport broadening by considering interface rough- 
ness (IFR) scattering with A = 4 A and A = 43 A, LO phonon scattering, and LA phonon scattering. The contribution 
of alloy disorder scattering was absent because the GaAs QW had AlAs barriers, and the influence of ionized impurity 
scattering was sufRciently reduced by the spacer layers. 

The calculated results for the lincwidth 2rop versus temperature are also shown in Fig. || by dashed (IFR), dash- 
dotted (IFR-fLO), and solid (IFR-f LO-I-LA) curves, in comparison with the transport broadening 2rtr. Additional 
width due to nonparabolicity is already included in these three curves, making small corrections compared with the 
contribution of interface roughness scattering, as seen in Fig. ^ At low temperatures, interface roughness scattering 
contributes 10.4 meV to lincwidth, and LO-phonon spontaneous emission contributes 0.7 meV. Phonon scattering 
processes become more active as temperature rises, and LO and LA phonon scattering contribute 1.8 and 0.7 meV, 
respectively, to lincwidth at room temperature. These calculated results are in good agreement with the experimental 
data shown by solid circles. Note that the increase in dashed line (IFR) with increasing temperature is due to 
the nonparabolicity effect; the contribution of interface roughness scattering itself slightly decreases with increasing 
temperature, as expected from the energy dependence in Fig. |^. 

For the transport broadening 2Ptr, interface roughness scattering makes a dominant contribution of 0.73 meV at 
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FIG. 8: Intersubband absorption spectrum observed at 4.5 K in a GaAs/AlAs single QW. The linewidth 2rop was 11.1 meV 
and the transport energy broadening 2rtr was 1.2 meV, which was calculated from the mobility fi of 2.9 x 10'' cm^/Vs. 




FIG. 9: Temperature dependence of the intersubband absorption linewidth 2rop and transport energy broadening 2rtr (or 
mobility fi). Circles show experimental values, and lines show numerical results calculated by considering interface roughness 
(IFR), LO phonon, and LA phonon scattering. 



low temperatures, which nearly explains the experimental value of 1.2 meV. As already pointed out in the previous 
section, this value of 0.73 meV is an order of magnitude smaller than the contribution of 10.4 meV to linewidth. In the 
temperature range above 80 K, the contribution of LO phonon scattering to transport broadening rapidly increases as 
temperature rises, and reaches a dominant value of 9.3 meV at 300 K, as is well known. Such an effect of LO phonon 
scattering on transport broadening is very different from that on linewidth. The contribution of LA phonon scattering 
to transport broadening is 1.2 meV, which is comparable with that to linewidth. 

As a result, linewidth and transport broadening have very different dependences on temperature. Similar behavior 
of linewidth versus temperature was also reported for GaAs/Alo.sGao.yAs QWs by AUmen et al.^. 



B. Well- width dependence 



Interface roughness scattering is expected to give absorption linewidth considerably strong well- width dependence. 
Campman et al. reported low-temperature linewidths and mobilities in modulation-doped GaAs/Alo.sGao.yAs QWs 
with Ng ~ 6 X 10 cm~ for various well widths in the range L = 75- 110 A g. Here, we calculate linewidth and 
transport broadening for the same structures. As scattering mechanisms, interface roughness scattering with A = 3 A 
and A = 85 A and LO phonon scattering are included one by one. 

Figure |l^ shows the calculated results for low-temperature linewidth versus well width in the range L = 75- 110 A. 
First, the well- width dependence of linewidth due to interface roughness scattering, shown by the dashed curve (IFR), 
is not so strong for small L, because the confinement of the first excited state is weaker and thus Fn = dEi/dL is 
considerably smaller in GaAs/Alo.sGao.yAs QWs than in infinite-barrier QWs. Second, the contribution of LO phonon 
scattering slowly increases as QWs become wider, which makes the well- width dependence of linewidth slightly weaker. 
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FIG. 10: Well- width dependence of intersubband absorption linewidth, calculated at OK by considering interface roughness 
(IFR) and LO phonon scattering. Solid circles show experimental results measured at low temperatures by Campman et al. 
[9]. 
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FIG. 11: Well-width dependence of transport energy broadening, calculated at OK by considering interface roughness (IFR) 
scattering. Open circles show experimental results measured at low temperatures by Campman et al. [9]. 



The solid curve (IFR+LO) is in good agreement with experimental results shown by solid circles If barriers are 
higher, as in GaAs/AlAs QWs, the first excited state is more strongly confined and the interface roughness scattering 
contributes much more to linewidth than LO phonon scattering does, which will lead to a much stronger well-width 
dependence of linewidth. 

On the other hand, the well-width dependence of low-temperature transport broadening is shown in Fig. |l^. The 
transport broadening considered here is determined only by interface roughness scattering, because intrasubband LO- 
phonon emission and absorption are impossible at low temperatures. Foo^ = (dEo/dL)'^ in Eq. ( |27| ) is proportional to 
in the infinite-barrier approximation, and this leads to a strong well-width dependence of transport broadening 
even in finite-barrier QWs. The calculated curve explains the experimental results plotted by open circles very 
well. 



C. Alloy composition dependence 

Experimental measurements of the alloy composition dependence of linewidth are expected to show the effects of 
alloy disorder scattering. Campman et al. reported low-temperature linewidths and mobilities in modulation-doped 
Ina;Gai_a;As/Al(, gGao.rAs QWs with L — 100 A and Ns 8 x 10^^ cm~^ for various compositions in the range 
X = — 0.1 H. We calculate linewidth and transport broadening for the same structures. As scattering mechanisms, 
interface roughness scattering with A = 3.5 A and A = 40 A, LO phonon scattering, and alloy disorder (AD) scattering 
are included one by one. 
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FIG. 12: Alloy composition dependence of intersubband absorption linewidth, calculated at OK by considering interface rough- 
ness (IFR), LO phonon, and alloy disorder (AD) scattering. Solid circles show experimental results measured at low tempera- 
tures by Campman et al. [9]. 
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FIG. 13: Alloy composition dependence of transport energy broadening, calculated at K by considering interface roughness 
(IFR) and alloy disorder (AD) scattering. Open circles show experimental results measured at low temperatures by Campman 
et al. [91. 



Figure shows the calculated results for low-temperature linewidth versus alloy composition x in the range 
X — — 0.1. The contribution of interface roughness scattering is 1.6 meV at a; = 0, and slowly increases as x 
increases because QWs become deeper. LO phonon scattering contributes approximately 1 meV to linewidth, almost 
independently of x. Although the contribution of alloy disorder scattering is proportional to x for small x, it is as small 
as 0.24 meV even at a; = 0.1. Our calculations explain the experimental observation of linewidth being insensitive to 
alloy composition, plotted by solid circles |9| . 



On the other hand, transport broadening is shown in Fig. 13 as a function of x. Interface roughness scattering 
contributes 0.1 meV to transport broadening, while alloy disorder scattering makes the larger contribution of 0.27meV 
at X — 0.1; this shows that transport mobility drops remarkably as x increases. Our calculations explain the exper- 
imental results plotted by open circles The small disagreement may be due to clustering in alloy layers, where 
the effective correlation length of alloy disorder in terms of roughness scattering may be larger than aj \p2/K in actual 
samples grown by molecular beam epitaxy (MBE) or metalorganic chemical vapor deposition (MOCVD). 

It should be noted that the one-order-of-magnitude different contributions of interface roughness scattering to 
linewidth and transport broadening are important in explaining their different behaviors versus alloy composition. 
Alloy disorder scattering in itself contributes fairly equally to linewidth and transport broadening, as mentioned in 
the previous section. 
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D. Doping position dependence 

If the donor doping position is varied, then the contribution of ionized impurity scattering to hnewidth should 
change. Dupont et al. measured low-temperature linewidths in J-doped GaAs/Alg 25Gao.75As QWs with L = 76 A 
and Ng ~ 1 x 10^^ cm~^ for two different doping positions: Z — and 1 12 A ||T^ . We calculate linewidth and transport 
broadening for the same structures. As scattering mechanisms, interface roughness scattering with A = 5.66 A (2 
MLs) and A = 70 A, LO phonon scattering, and ionized impurity (ION) scattering are included one by one. 

Figure |lj shows the calculated results for low-temperature linewidth versus doping position Z in the range Z = — 
120 A. Interface roughness scattering and LO phonon scattering contribute 5.8 and 0.8 meV to linewidth, respectively. 
When donors are doped in barriers, at Z ^ 100 A for example, the contribution of ionized impurity scattering is as 
small as 0.3 meV. When donors are doped in QWs, at the center Z = oA for instance, the contribution of ionized 
impurity scattering is 2.8 meV, which is smaller than that of interface roughness scattering. Our calculations explain 
the experimental results plotted by solid circles [0. 

Note that the wave function Ci(^) of the first excited state penetrates largely into the low barriers in these narrow 
GaAs/Alg gsGao.TsAs QWs, so the effect of ionized impurity scattering is greatly enhanced even in barrier-doped 
QWs. If wave functions are more strongly confined, for example, as in the narrow GaAs/AlAs QWs used in our 
experiment, the contribution of ionized impurity scattering to linewidth is less than 0.1 meV for barrier-doping. 

On the other hand, low-temperature transport broadening is shown in Fig. |l^ as a function of Z. Interface roughness 
scattering contributes 0.44 meV to transport broadening, while ionized impurity scattering contributes 12.2 meV at 
Z — A. and 0.33 meV at Z — 100 A. Therefore, mobility greatly decreases when donors are doped in or near QWs; 
rather thick spacer layers, more than 150 A in this case, are necessary to completely remove the influence of ionized 
impurity scattering on mobility. 

IV. SUMMARY 

We have formulated the microscopic energy-dependent relaxation rate 2rop(-E) of intersubband optical transition 
in QWs due to scattering by interface roughness, LO phonons, LA phonons, alloy disorder, and ionized impurities, 
and have numerically calculated the absorption linewidth 2rop for GaAs-based QWs in comparison with the transport 
energy broadening 2rtr = 2he/m*fj, related to the mobility /i. 

The sensitivity of linewidth to interface roughness scattering is about one order of magnitude higher than that 
of transport broadening, because the contribution from the intrasubband scattering in the first excited subband is 
larger than that in the ground subband. This provides an essential insight for understanding experimental values for 
linewidth and the apparent lack of correlation between linewidth and mobility. 

The contribution of LO phonon scattering to linewidth is small, about 2 meV in narrow GaAs-based QWs even at 
room temperature, because the difference in intrasubband scattering matrix elements for the two subbands is small 
owing to the cancellation of form factors. In addition, intersubband LO-phonon spontaneous emission contributes 
approximately 1 meV to linewidth at low temperatures. Therefore, linewidth has a very weak temperature dependence, 
while mobility is greatly lowered by LO phonon scattering in the temperature range above 80 K. 

LA phonon scattering and alloy disorder scattering give matrix elements that are independent of scattering vectors, 
and lead to a linewidth comparable with transport broadening. The contribution of LA phonon scattering is, for 
example, about 1 meV at room temperature in narrow GaAs-based QWs, and this is small for linewidth compared 
with the contribution of interface roughness scattering. Alloy disorder scattering contributes, for instance, about 
0.3 meV in In^^Gai-ajAs QWs with x — 0.1. This is negligible for linewidth but predominant for transport broadening, 
causing a remarkable drop in mobility as x increases. 

Ionized impurity scattering contributes little to linewidth in modulation-doped QWs, because the difference in 
intrasubband scattering matrix elements for the two subbands is small. On the other hand, rather thick spacer layers, 
more than 150 A in narrow GaAs/Alo.sGao.rAs QWs for example, are necessary to remove the influence of ionized 
impurity scattering on mobility. 
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FIG. 14: Doping position dependence of intersubband absorption linewidth, calculated at K by considering interface roughness 
(IFR), LO phonon, and ionized impurity (ION) scattering. Solid circles show experimental results measured at low temperatures 
by Dupont et al. [10]. 
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FIG. 15: Doping position dependence of transport energy broadening, calculated at OK by considering interface roughness 
(IFR) and ionized impurity (ION) scattering. 



TABLE I: Material constants of GaAs. 



band gap of Al^^Gai-x As (x < 0.45) at K (1.519 + 1.247a;) eV 

band gap of AlAs at K 3. 1 13 eV 

conduction-band discontinuity ratio for GaAs/AlGaAs ~ 0.65 

static dielectric constant kq = 12.91 

optical dielectric constant Koo = 10.92 

LO phonon energy ?iwlo = 36.5 meV 

deformation potential constant D — 13.5 eV 

longitudinal elEistic constant c; = 1.44 x 10^^ N/m^ 

spin-orbit splitting 0.341 eV 

Kane energy 22.7 eV 
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